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ABSTRACT 


Results of an experimental study on the dynamic properties of 
tantalum are presented. Areas studied included stress-strain- 
strain rate and reverse loading behavior, elastic constants, 
equation of state, compressive and release wave characteristics, 
and spall fracture. The material showed approximately elastic- 
plastic behavior with some strain hardening under uniaxial 
stress compression, and exhibited strain-rate sensitivity with 
yield under uniaxial stress increasing from 1.4 to ^ 5.4 kbar 
for a strain rate increase of 0.001 to 800/sec. Longitudinal 
and shear wave velocities at 20°C were 4.146 and 2.032 mm/ysec, 
and temperature and pressure dependence was also measured and 
various elastic constants were calculated. The shock wave equa¬ 
tion of state (EOS) up to 200 kbar was determined as a H = 4.0 + 

550 u + 215 u 2 . The EOS was also calculated from the elastic 
P P 

constants and comparisons made with the shock wave EOS. Com¬ 
pressive wave tests showed a well-defined elastic precursor with 
a steady-state value of 'v 15 kbar. The impact velocity required 
for spall fracture was found to increase v/ith decreasing pulse 
width. 
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FOREWORD 


The prediction of reentry vehicle response to impulsive loading 
resulting from energy deposition has been studied extensively 
during the last decade. Analytical models and computer routines 
have been developed to assess the vulnerability of such vehicles 
to the initial loading phase, in which stress waves are gener¬ 
ated and propagated through the structure, and to subsequent 
phases, where elastic vibration, plastic deformation and/or frac¬ 
ture may occur. 

A Defense Nuclear Agency (DNA) program was initiated with the ob¬ 
jective of improving computer codes for the prediction of damage 
induced by X-rays (PREDIX)» The PREDIX metals program has combined 
the efforts of several contractors under the direction of DNA. 

The primary contributions of each contractor can be summarized as 
follows: 

General Motors Corporation—Measurement of material 
properties, including stress-strain-strain rate be¬ 
havior, elastic constants, equations of state, shock 
wave profiles and spall fracture. 

Effects technology Inc.—Spall fracture tests and 
analysis (exploding foil) and underground testing. 

KMS Technology Center—Constitutive modeling. 

Physics International Company—Energy deposition and 
spall fracture studies (electron beam). 

Systems, Science and Software—Energy deposition and 
front surface phenomena, constitutive modeling and 
code development, and underground testing. 
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INTRODUCTION 


This report is one of a series of six giving results of a study 
of dynamic properties of four metals. The principle objectives 
of this study were: 

• Provide measurements of dynamic properties of 
materials to be used as inputs to model develop¬ 
ment for shock wave propagation and dynamic 
fracture. 


• Provide data to establish accuracy of computer 
code predictions. 

• Provide physical interpretation of experimentally 
observed material response to dynamic test condi¬ 
tions . 

The primary areas studied included: (1) Stress-strain-strain 
rate behavior in uniaxial stress, including Bauschinger effect; 
(2) Elastic constants, including temperature and pressure 
dependence; (3) Hugoniot equation of state; (4) Compressive 
and release wave profiles, including elastic precursor decay 
and wave attenuation; and (5) Spall fracture studies, including 
pulse width dependence. The six reports in this series are: 


Vol. 

I: 

Summary of Results, DASA 2501-1 (AD712847) 

Vol. 

II: 

Experimental Methods and Techniques, 



DASA 2501-2 (AD 730750) 

Vol. 

III: 

6061-T6 Aluminum, DASA 2501-3 

Vol. 

IV: 

Alpha Titanium, DASA 2501-4 

Vol. 

V: 

OFHC Copper, DASA 2501-5 (AD 728846) 

Vol. 

VI: 

Tantalum, DASA 2501-6 
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The experimental techniques utilized in this study are reviewed 
in Vol. II (DASA 2501-2). Briefly, the principal techniques 
applied were: 

Stress-strain Studies—A laboratory-type universal 
testing machine was used for rates < 0.1/sec, a 

_3 

medium strain pate machine for rates of 10 /sec 
2 

to 10 /sec, and a Hopkinson bar device for rates 
of about 10 2 to 5 X 10 2 /sec. 

Elastic Constants—The pulse superposition method was 
used for measuring longitudinal and shear wave velo¬ 
cities as functions of temperature and hydrostatic 
pressure. 

Equation of State—Gun-launched, flat-plate impact tech¬ 
niques were used to generate uniaxial strain conditions. 
Velocities up to 0.6 mm/ysec were achieved with 63.5 mm 
and 102 mm single-stage compressed-gas guns. Hugoniot 
data were obtained with x-cut quartz gages. 

Wave Prot agation and Spall Fracture—Flat-plate impact 
techniques wer^ also used. Wave profiles were measured 
with x-cut quartz gages and a laser velocity interfero¬ 
meter. 

A literature search on dynamic properties of materials was 
carried out under the PREDIX program and a listing of recent 
publications {including abstracts) on dynamic properties of 
materials is given in Reference 1. 


2 
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SECTION I 

MATERIAL PROPERTIES 


Tantalum has a body-centered-cubic structure and the material 
on which tests discussed in this report were carried out was 
a commercially pure £i^de of 99.5% purity. This material was 
purchased as 50.8 mm diameter bar stock, which was preferrable 
to plate or sheet stock since it provided specimens of various 
thicknesses and gave little scrap while ensuring that all speci¬ 
mens were from the same parent material. Chemical analysis 
performed by the supplier* gave: 


Columbium 

210 

Tungsten 

1 

Iron 

20 

Carbon 

<10 

Oxygen 

38 

Nitrogen 

22 

Hydrogen 

<5 


Tantalum is not heat-treatable and different strength levels 
are achieved by slight changes in impurity content or by cold¬ 
working. Although not specified, the material apparently had 
been annealed. Certifications gave 1.8 kbar tensile yield 


*Fansteel, Inc., Lot No. 69B-2522, 31 July 1969. 
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strength, 2.35 kbar ultimate tensile strength and 52% elonga¬ 
tion at 4 d (measured along the bar axis). Measured hardness 
was 52 A photomicrograph of a representative portion of 

the bar is shown in Figure 1. The bar axis or centerline is 
vertical and the width shown is ^ 5 mm out of a total o^ 

50.8 mm. There was a slight grain elongation along the axis 
and about a factor of 10 variation in grain size, ranging from 
ASTM 6 to 00 (45 to 500 vim). For results presented in this 
report, the test or wave propagation direction was axial 
(e.g., vertical in Figure 1). 

3 

The average measured density (p Q ) was 16.66 g/cc (0.601 lb/in ). 

For use in equation of state and energy deposition claculations, 

(2-5) 

several physical constants were compiled from the literature 
and are listed below: 

■ 

Volume coefficient of thermal expansion, B: 

19.8 x 10" 6 + 0.001 x 10" 6 T/°C 
Specific heat, c^: 

0.034 + 10 -5 T cal/g/°C 
Melting temperature: 2996°C 
Boiling temperature: 5430°C 
Latent heat of fusion: 35 cal/g 
Latent heat of vaporization: 1000 cal/g 
(T in °C) 

A portion of the tantalum bar was annealed for 1 hour at 1200°C 
in vacuum. The structure of this material is compared to the 
as-received material in Figure 2. There was no significant 
change in structure or hardness. 


4 
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SECTION II 

STRESS-STRAIN STUDIES 


Shock wave and structural response code input requirements 
include constitutive equations relating stress to strain and 
strain rate. Uniaxial stress tests at various strain rates 
give yield and flow stress behavior, strain rate sensitivity 
and work hardening characteristics. Bauschinger effect tests 
were conducted to provide data on unloading and subsequent 
yield behavior. 


UNIAXIAL STRESS TESTS* 


The compressive stress-strain-rate behavior of tantalum at 20°C 
is shown in Figure 3 for strain rates of 10 /sec to 600/sec, 
with each curve the average of three tests. The material is 
very strain-rate sensitive, showing an almost 300% increase in 
yield level over this range of strain rates. A slight yield 
drop is evident above a rate of 40/sec. Strain rate data are 
cross-plotted in Figure 4 to give true stress vs. log true 
strain rate. Both yield and flow stress (6% strain) increase 
with strain rate, with sensitivity (slope) becoming very high 
above 100/sec. 


* Deformation mechanisms of tantalum (primarily under uniaxial 
stress conditions) have been studied by a number of other 
investigators. The reader is referred to the literature for 
details in such areas as strain rate effects (Refs. 6-10, 
temperature effects (0-14), impurity content (6, 13-16), 
stress relaxation (9,17) grain size effects (14, 17) disloca¬ 
tion behavior (11, 15-18), anisotropy (IS), and fracture (12). 


6 
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Figure 3 Compressive Stress vs. Strain, Tantalum 



Figure 4 Compressive Stress vs. Log Strain Rate, Tantalum 
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REVERSE LOADING TESTS 

The reverse loading behavior or Bauschinger effect was studied 
by performing uniaxial tension tests after the material hac been 
prestrained in the opposite (compressive) direction. Results 
from low strain rate (0.001/sec) tests are given in Figure 5 
for maximum strains of 0.5, 1, and 3%. The Bauschinger strain 
(defined here as the plastic strain in the reverse loading path 
at 3/4 of the initial yield stress) is insensitive to prestrain, 
at least up to 3% strain. Flow stress in tension has increased 
with increasing compressive prestrain. 



Figure 5 Reverse Loading Behavior, Tantalum 
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Cyclic loading*behavior is shown in Figure 6 where a single speci 
men was put in compression and then tension through a total of 
5 cycles. The observed yield and flow stress behavior is indica¬ 
tive of an isotropic hardening mechanism. 


M j I • — 

-1.4 


J. 

.4 



Figure 6 


Cyclic Loading Behavior, Tantalum 
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SECTION III 

ELASTIC BEHAVIOR 


Measurement of elastic constants provides additional inputs to 
the study of material response. Three basic measurements are 
density, longitudinal (dilatational) wave velocity and shear 
(transverse) wave velocity. Measurements of temperature and 
pressure dependence of the wave velocities lead to prediction 
of pressure-compression isotherms, adiabats and hydrostats. 

If these measurements are made with sufficient accuracy (%0.1%), 
the calculated pressure-compression behavior can frequently be 
extrapolated with reasonable confidence to higher pressures 
than covered in the actual measurements. 

WAVE VELOCITY MEASUREMENTS 

Initial measurements of longitudinal and shear wave velocities 
were made at 20°C and atmospheric pressure (i.e., P = 0) . 

Values obtained parallel to the bar axis were: 

C_ = 4.146 + 0.005 mm/ysec 

Xj — 

Cg = 2.032 + 0.003 mm/ysec 

Although the accuracy of these measurements are estimated to be 
about +0,1%, it should be noted that C L and C g may vary depending 
on the structural form of the material (i.e., plate, bar, etc.), 
its thermal and mechanical history, and the measurement direc¬ 
tion relative to any structural anisotropy. 
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Measurements were made of temperature and hydrostatic pressure 
dependence of wave velocities and the data corrected for thermal 
expansion and hydrostatic .compression *to give: 

\ . I 

1 ‘ ! 

Temperature Dependence, 0-100°C, P = O' 

, I 

C T = 4.149 - 0.000146T mm/ysec 

1 

C c = 2.036 - 0.000173T mm/ysjec 

b * , I 

(T in °C) , 1 i 

i 

i ; 

Pressure Dependence, 0-9 ‘kbar, T = 20°C 

C T = 4.146 + 0.00291P mm/usec . , , 

Jj * i 

C c = 2.032 + 0.00137P mm/ysec 

' I » 

(P in kbar) , 1 

: , . 1 

Temperature and pressure dependence of elastic 'constants for 

( 20 ) 1 

single-crystal tantalum have been reported by Soga and. 

pi \ (22) ( 2 3 i 

Chechile, * and Palmieri' ' and Armstrong. and ,BrQ’:/n v ' , 

i 

have published data on temperature dependence of elastic con- 

, , i 

stants in polycrystalline tantalum. Palmieri also found C„ 

f % & S 

(sound wave velocity) to vary from 3.245 to 3.366 mm/ysec and, 

C_ from 2.035 to 2.085 rom/ysec (at 0°C), depending on specimen 

(24) 

orientation and heat treatment (i.e., annealing)) Lamberson 
studied both temperature and pressure dependence of C L and C g 
in polycrystalline tantalum with a density of 16.64 g/cc and 

I 

16 to 30 ym grain size. was given as 4.156-0.0001573T and 
4.153 + 2.0025?, and C g as 2.037 - 0.0001418T and 2.034 +, 0.0 ! 010P 
(T in °C and P in kbar). ‘ 


11 
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. . i 

ELASTIC CONSTANTS 

* ’ ' 

.Assuming polycrystalline tantalum is isotropic (C T and C 

' | iJ u 

are independent 6f direction of measurement in a polycrystal¬ 
line specimen), wave velocity data can be used to calculate 

■ : i 

various elastic constants. At 20°C and zero pressure, the 

i ‘ 

following adiabatic constants.were obtained: 

! i : 

! ! J 

Bulk wave velocity, C fi = 3.42'mm/ysec 
Sound wave velocity, C_ = 3.33 mm/ysec 

« , £j 

i Rayleigh wave*velocity, C D = 1.90 mm/ysec 

I j f K 

Poispon's ratio, v'=.0.3&2 
1 Bulk modulus, K = 1946 kbar 
•. Shear modulus,- G = 688 kbar 

' 1 .i 

■ Elastic modulus, E - 1846 kbar 

• » * r 

Lame's parameter, A = 1488 kbar 


I 


I 


• * 

The adiabatip bulk and shear, moduli were obtained as functions 
of T and P from the following: 



G J pc/ 


I 


* A cubic material is 'elastically isotropic when — = 1, 

C 11“ L 12 

where |C 4< |, C^ and C^ 2 are elastic stiffness constants. Single- 
= crystal tantaluirij is slightly anisotropic with — = 

' or orv C ll“ C 12 

1 >56 (20,25,26). 
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This gave: 

Temperature dependence, P = 0 
K S = 1948 - Q.C84T kbar 
G S = 691 - 0.131T kbar 

Pressure dependence, T = 20°C 

K S = 1946 + 3.79P kbar 
S 

G = 688 + 1.28P kbar 

where the superscript S indicates adiabatic. 

Isothermal values of K and G as well as adiabatic and isothermal 
pressure and temperature derivatives at 20°C and zero pressure 
were calculated and results are given in Table I. Use of these 
constants in calculating the isotherm, adiabat and hydrostat is 
discussed in the section on Equation of State. 


GRUNEISEN PARAMETER 

The Gruneisen ratio y is a parameter in the solid equation of 
state relating pressure to volume and energy, y = V (3P/3E) V 
(see, e.g.. Ref. 27). This parameter can be expressed thermo¬ 
dynamically as: 



13 
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TABLE I 

ELASTIC CONSTANTS FOR TANTALUM 
(P = 0, T = 20°C) 


PARAMETER 

VALUE 

K S 

1946 kbar 

= K S’ 

V 3P 'T 0T 

3.79 

( H _) p 

- 0.084 kbar/°C 

<#> s - * 

3.77 

( 55 ®, 

'3T 1 S 

15.1 kbar/°C 

T 

K 1 

1928 kbar 

tif.) = k t ' 

'3P ; T 0T 

3.80 

& 

V 3T 'P 

- 0.146 kbar/°C 

(*£.) = k t ’ 

'3P ; S ^OS 

3.77 

T 

/3R v 
{ W> s 

15.1 kbar/°C 

G S - G T 

688 kbar 

,££. _ ' 

'3P'T 0T 

1.28 

( 3 £, 

V 3T'P 

- 0.131 kbar/°C 

_ • 

'3P'S “ ''OS 

1.25 


5.02 kbar/°C 
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For an isotropic elastic solid, this gives: 



where C T and C.. are measured under adiabatic 

Jj O 

conditions. 

At 20°C and zero pressure; 

Y 0 - 1.63 

The zero-pressure Gruneisen parameter can also be estimated by 

(28) 

yeral other methods, including those of Slater , Dugdale 

(29) (30 31) 

c»nd MacDonald , Anderson and Dienes ' , and Schuele and 

(30 32) 

Smith ' . These methods are. discussed in Appendix A. 

An estimate of the temperature dependence of the Gruneisen 
parameter can be obtained by differentiating Equation 1 with 
respect to temperature. At constant pressure (P = 0, T = 20°C): 



The pressure dependence can be estimated by differentiating 
with respect to pressure. At constant temperature (T = 20°C, 
P = 0): 
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Another approach is to assume y/V constant, which permits 
Oy/3P)„ to be estimated directly from the bulk modulus data. 
At constant temperature (20°C): 



= -11.93 

(1928+3.8P) 1,263 

= -0.00085/kbar @ P = 0 


DEBYE TEMPERATURE 

The Debye temperature 0 is important in thermal energy calcu¬ 
lations and indicates the temperature above which variations 
due to temperature for some thermodynamic parameters such as 
specific heat, thermal expansion and Gruneisen ratio become 

small. There are a number of methods for calculating the 

(33) 

Debye temperature , and one approximation suitable for use 
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with elastic wave velocity data for polycrystalline metals is: 



where h is Planck's constant, k is Boltzmann's 
constant, N is number of mass points and V is 
sample volume. 

At 20°C, the elastically-determined Debye temperature for 
tantalum is: 



0 = 260°K 


-7 3 

for N = 2 (atoms per unit cell), V = (3.306 x 10 mm) 

(34) 

(lattice constant) , C L = 4.146 mm/usec and C g = 2.032 
mm/usec. 
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SECTION IV 

EQUATION OF STATE 


The development of material response models and computational 
codes requires data on material behavior under uniaxial strain 
conditions/ including equation of state, wave propagation and 
spall fracture. Most data are used directly in developing 
models of material behavior, although independent check data 
(attenuated wave and spall wave profiles) are necessary to 
determine accuracy of the calculations. This section presents 
equation of state analyses, and the following two sections 
cover wave propagation and spall fracture. 


SHOCK WAVE EQUATION OF STATE 

The hugoniot equation of state is the locus of equilibrium 
states reached after shocking of a material, and data are 
usually obtained either as stress-particle velocity points or 
as shock velocity-particle velocity points. The hugoniot data 
presented in this report were obtained with x-cut quartz 
gages and representative records for direct impact and trans¬ 
mitted wave tests are shown in Figure 7. The buffered direct 
impact method (tungsten carbide buffer plate on front of the 
quartz) permitted stresses up to 80 kbar in tantalum while 
keeping the stress in quartz at an acceptable level. The use 
of quartz gages for transmitted wave tests was primarily for 
the study of compressive wave development and elastic precursor 
decay, and results are discussed in the Wave Propagation section. 
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DIRECT IMPACT 
Ta-*Q 

Vj» 0.164 mm/ /is 

O'" 20. 8 kbar 
H 



DIRECT IMPACT, BUFFERED 
Ta-^WC/Q 
V,» 0.159 mm//is 
O'- 60.8 kbar 



TRANSMITTED WAVE 
AI-^Ta/Q 
V |- 0.269 mm/ /is 
O’" 32.8 kbar 


Figure 7 Quartz Gage Records 
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The direct impact records showed very fast rise times (<10 nsec) 

and then a straight peak stress region. The slight rounding at 

the front of the direct impact records in Figure 7 is due to 

impact tilt, and the ramped peak stress region results from 

finite-strain effects in the quartz which are corrected in the 
(35) 

data analysis. 

The hugoniot may be expressed in several forms. A convenient 
form for experimental work is that established by a least-squares 
fit to data in the stress-particle velocity (a H - u^) plane. 
Transformation of the hugoniot into other planes, such as shock 
velocity-particle velocity (U g - Up) or stress-volume (a H - v), 
is performed by assuming a material model. This was done by 
assuming an ideal elastic-plastic wave structure with equilib¬ 
rium initial and final states and applying the mass and momentum 
conservation equations: 


" °e + °e (D S - u e> (u p ' u e> 


and 


v = 


v. 




V»e ) 


(7) 


( 8 ) 


where p g and u g are density and particle 
velocity at the elastic limit c e . 

Therefore, when either a H - u p or U g - u p relations are estab¬ 
lished and a , p . u and C T are known, a„, U_, u and v are 
e e e L H s p 

uniquely determined. The resulting hugoniots for tantalum at 
20°C are given in Figures 8, 9, and 10, and are listed below: 
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0 0.10 0.20 0.30 0.40 


PARTICLE VELOCITY (mm/ us) 

Stress-Particle Velocity Hugoniot, Tantalum 



0 0 . 02 0.04 0.06 0 . 08 0 . 10 


COMPRESSIONS' hj 2 - ' l 


Figure 9 


Stress or Mean Pressure-Compression 
Hugoniot, Tantalum 
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Figure 10 Shock Velocity-Particle Velocity 
Hugoniot, Tantalum 


o H = 4.0 + 550 Up + 215 u 2 (Std. Error 0.3 kbar) 
o„ = 7.2 + 1820 ji + 3120 u 2 

n 

P fi = 1820 v> + 3120 u 2 (Hydrostat) 

U e = 3.36 + 1.239 u 
S p 



J 


where p = (v /v - 1) 
o 

These equations are based on the data points indicated by 

the symbols in Figure 8 and on a = 20 kbar/ and are eon- 

e * 

sidered to be valid over the stress range of 20 to 200 kbar. 
Accuracy of the quartz gage data is estimated as +2% for 
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o„ and + 1% for u . Hugoniot points reported by Rohde and 
n (27) p 

Towne for 99.9% polycrystalline tantalum at 30.2, 43.1, 

48.8 and 85.8 kbar are within +1%, - 8% of the fit to the 
data given in Figure 8. 

The hydrostat ( p H “V0 was obtained by subtracting the stress- 
axis intercept value from the ct h ~y equation. This approach 
satisfies- the requirement that the P^-p function pass through 
the origin and assumes the deviatoric stress (a D ) is constant 
(no strain hardening). Alternative ways of estimating 
for the low-pressure region are discussed in Appendix B. 


ULTRASONIC EQUATION OF STATE 


The equation of state can also be determined from ultrasonics 
measurements. Accurate measurement of elastic wave velocities 
(and, therefore, bulk modulus) as a function of hydrostatic 
pressure makes it possible to directly estimate shock wave com¬ 
pression behavior. Numerous analytical and empirical relations 
have been developed for relating pressure, volume and bulk 
modulus, including those of Birch (38*39), jjurnaghan^^ and 
Keane^^'^. The use of these relations is discussed in 
Appendix C and results are summarized below: 


Birch EOS 








P T = 2891 

U + m> 7/3 - d + y) 5/3 


1-0.15 

(1+y) 2/3 -1 



- 


_ 




The hugoniot was extended to 200 kbars by including three 
calculated points obtained from the high-pressure U_-u 
hugoniot given in Reference 36. This hugoniot was converted 
to the plane after assuming a hugoniot elastic limit 

of 20 kbar, and three 0 H “Ug "points" were then calculated 
at particle velocities of U.20, 0.25 and 0.30 mm/psec. 
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Murnaghan EOS, Isotherm 


P T = 507.1 


(1+y) 3,80 -1 


Murnaghan EOS, Adiabat 


P s = 516.2 


(i+y) “1 


Keane EOS, K S ' = 3.6 
00 


P s = 566.1 


(1+y) 3 * 6 -1 


- 91.9 In (1+y) 


( 11 ) 


( 12 ) 


(13) 


Keane EOS, =1.10 


Pg = 6063 


(1+y) 1 * 10 - 1 


- 4723 In (1+y) 


(14) 


In the above equations, subscripts T and S indicate isothermal 
and isentropic, respectively. 


Duvall gives a method of calculating the increase in entropy 

across a shock which can be applied to the Murnaghan equation 

. . (24,43) 

to give: 


Y 0 K 0 

P = P 4- —__ 

*11 *S 12 


P H = P s + 1261 


( K 0S + X )(^) 

(*)’ 


(15) 
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; 

Zel'dovich gives an expression relating P g and assuming 
Y/v is constant: ^ 



The above equations of state are Compared to the -shock wave 
EOS in Table II. The entropy correction based on the JMufnaghan 
equation has been applied to the Birch and Keane equations as. 
well. The ultrasonic equations of ; state show good agreement 
with each other but are 3 to 6% above the shock wave hydrostat. 

• i 

■ 

The shock velocity-particle velocity relationship can also ( be 
determined from ultrasonic data. Following the method of 

(45) . ! ’ s 

Ruoff one obtains: 



i 


i 
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For tantalum: 

U_ = 3.42 + 1.193 u + 0.14 u 2 
S P P 

2 S" 

The u coefficient (A = 0.14) was calculated assuming K_ c = 0. 

P c i» Ub 

If one assumed A = 0, then IC c = -0.003. Since the uncertainty 
in K“ s is on the order of + 0.01, the A coefficient will be 
taken as zero, which gives: 

U_ = 3.42 + 1.193 u (19) 

b p 

This is compared with the linear U,, - u relation obtained 

S p 

from low-pressure shock wave data in Figure 10. 

Vaidya and Kennedy measured the hydrostatic compressibility 
of 99.99% polycrystalline tantalum rod up to 45 kbar.Re¬ 
sults were given as 

” 77 = l+F = 0-«53xl3‘ 3 P T - 0.59394 x10~ 6 P 3 (20) 

Comparing this with «.*e shock wave and ultrasonics results 
presented in this report gives: 

P T , ISOTHERMAL PRESSURE (kbar) 



Shock Wave 

Ultrasonics 

hydrostatic 

u 

Eg. 9* 

Eq. 11 

Eq, 20 

0.005 

9.2 

9.7 

10.2 

0.010 

18.5 

19.5 

20.5 

0.015 

27.9 

29.5 

31.0 

0.020 

37.5 

39.6 

41.6 


* 

Corrected for isothermal conditions. 


The hydrostatic data of Vaidya and Kennedy are 5 to 10% above 
(P^ - y plane) the shock wave and ultrasonic data. 
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YIELD BEHAVIOR 

A complete description of the low-pressure equation of state 
requires consideration of yield behavior or the hugoniot elas¬ 
tic limit. Values of compressive yield were obtained by three 
independent methods. 

A. Uniaxial Stress 


The yield level was determined in uniaxial stress 
(a^) as a function of strain rate and then converted 
to unaxial strain using: 


o 


e 



( 21 ) 


Strain Rate 
(sec" 1 ) 


800 

10" 3 


B. Wave Profiles 


(kbar) 


11.2 


2.9 


The elastic precursor level was measured as a function 
of propagation distance using quartz gages. These 
profiles are discussed in the section on Wave Pro¬ 
pagation. 


Propagation 

Distance °e 

(mm) (kbar) 


2 19 

12 16 
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C. Hugoniots 

The yield level in uniaxial strain was inferred by com¬ 
parison of the elastic (p rt C T ) and plastic (a„ vs. u ) 

U n p 

hugoniots, i.e., by defining the elastic limit as the 
intersection of the elastic response line and the curve- 
fit to the stress-particle velocity data 

10 P n C T u = A + Bu + Cu ^ 

0 L e e e 

Solving for u^ gives 

u = 0.029 mm/ysec 
e 

= 20.0 kbar 

For an isotropic, polycrystai'ine metal with elastic-perfectly 
plastic behavior and no strain-rate or time-dependent effects, 
one would expect the same yield to be obtained by each method. 
The yield obtained from wave profile tests and the hugoniot 
measurements are in fair agreement, but are both much higher 
than the yield determined from uniaxial stress data for rates 

3 

less than 10 /sec. For this particular tantalum, the equili¬ 
brium yield in uniaxial strain was assumed to be 20 kbar for 
use in the EOS transformations made in obtaining the a„ - y 

n 

and U c - u relations given above (Equation 9). 

p 
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t 

SECTION V 

WAVE PROPAGATION 


Profiles of shock waves propagated through a specimen were 

recorded as stress-time or velocity-time histories. The 

stress-time data were obtained with quartz gages and the 

results were transformed to material stress by application 

(35) 

of an impedance matching technique, assuming time-in¬ 

dependent behavior. The velocity-time data were obtained 
with a velocity interferometer and are presented as measured. 
All wave profiles obtained are given in Appendix D for ref¬ 
erence. 


COMPRESSIVE WAVE BEHAVIOR 

Structure in the compressive wave is shown in the quartz gage 

* 

data in Figure 11. The wave front is characterized by an 
elastic portion and a transition to a spreading plastic wave 
(nominal final stress of 34 kbar). (The elastic wave would 
be overdriven by the plastic wave at 'v 440 kbar) . 

Also shown in Figure 11 are shock velocities as calculated 
from the shock wave hugoniot (Equation 9). Although the cal¬ 
culated plastic wave velocity lies within the rise-time of 
the measured wave, the plastic wave is not a step pulse and 
shows spreading with propagation distance at this stress 


In this and subsequent figures containing quartz gage data, 
tilt refers to the time required for a step-input to sweep 
across the gage electrode diameter. 
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level. This means that transformation of the hugoniot from 

the o„ - u plane to a„ - u or U_ - u by assuming ideal elastic- 
H p H bp 

plastic wave structure and steady-state conditions may lead to 
some uncertainty at lower stresses. 



Figure 11 Compressive Waves, Propagation Distance 
Dependence (6061-T6 A1 Impactors) 

Two tests with quartz gages were conducted using specimens 

machined from a piece of 76.2 mm diameter bar, and results are 

compared to 50.8 mm bar data in Figures 12 and 13. Although 

the 76.2 mm bar was ordered to the same nominal specifications 

* 

as the 50.8 mm bar , the compressive wave development is signi¬ 
ficantly different. The larger diameter bar showed a less well- 
defined elastic limit but a "faster" plastic wave. These dif¬ 
ferences can probably be attributed to changes in metallurgical 
characteristics of the material. 


•/( 

Impurity content for the 76.2 mm bar was about the same as 
the 50.8 mm, but the yield and flow stress at low strain rate 
was about twice as high and hardness was 86 R„. 
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i 

I 



Figure 12 Compressive Waves, Material Differences 
(6061-T6 Al Impactors) 



Figure 13 


Compressive Waves, Material Differences 
(6061-T6 Al Impactors) 
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Precursor decay data for tantalum, as determined from quartz 
gage data, is summarized in Figure 14. The calculated initial 
elastic impact stress was 37 kbar for these tests and the 
precursor level has dropped ^ 50% in 2 mm of travel. The 
equilibrium or steady-state elastic limit appears to be *v» 15 kbar. 
The uncertainty bars in Figure 14 reflect dispersion and rounding 
at the elastic front which is due to at least three factors. 

First, since quartz gages average stress, over the electrode 
area, small differences in wave front arrival times at the 
specimen/gage interface would give a finite rise-time rather 
than an instantaneous stress jump. Second, the finite thick¬ 
ness of the epoxy between the specimen and the gage will in¬ 
crease the apparent rise-time in the wave front. Third, the 
influence of shock wave tilt on a finite-area gage is to 
smooth out abrupt changes in stress level as well as to in¬ 
crease recording rise-time of the wave front. 



Figure 14 Elastic Precursor Decay in Tantalum 
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Gillis, et al., ' reported elastic precursor decay data for 
annealed, high-purity (99.9%) tantalum. Direct comparison with 
the present work is not possible because of material differences 
and variations in impact stress and propagation distance, 
however, their data tend toward the same steady-state yield 
value, 15 kbar. Gillis, et al., used quartz gages to record 
the stress profits and indicated elastic wave rise-times on 
the order of 100 nsec, after correction for tilt. The quartz 
gage data shown in Figure 11 ( not corrected for tilt), show 
rise-times generally less than 50 nsec, and velocity inter¬ 
ferometer data, discussed below, showed rise-times to be less 
than 20 nsec. 

To determine if rise-time of the input pulse influenced pre¬ 
cursor development, a series of velocity interferometer tests 
were performed with a fused quartz buffer on the front of the 
tantalum target. (Test conditions for all velocity inter¬ 
ferometer tests discussed in this report are listed in Table III). 
Because of the nature of the fused quartz hugoniot at stresses 
below 'v 40 kbar, the compressive wave spreads as it 
propagates. The rise-time of the pulse at the fused 
quartz/tantalum interface was ^300 nsec and the resulting 
compressive wave profiles after propagation through several 
thicknesses of tantalum are shown in Figure 15. The calculated 
elastic limits at 6 and 10 mm (18 and 16 kbar, respectively) 
are aDoroximately the same as those determined from the quartz 
gage tests. This indicates that elastic wave development is 
insensitive to input pulse rise-time (i.e., strain rate), at 
least for rise-times of 300 nsec or less. Note that Test 225 
in Figure 15 shows evidence of a slight yield drop or stress 
relaxation. Although this behavior might "be expected on the 
basis of the relaxation observed in the uniaxial stress tests 
(see Figure 3), it was seen on only two velocity interferometer 
tests and on none of the quartz gage tests. 
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^'■^•Uure 15 . Compressive Waves , Long Rise-Time 

Input Pulse 


RELEASE WAVES AND WAVE ATTENUATION 

Release waves were studied by using relatively thin impactors 
and measuring the complete wave profile with the laser velocity 
interferometer. Unattenuated wave profiles are shown in 
Figure 16 for unbacked (free rear surface) and Lexan-backed im¬ 
pactors. The unloading wave in tantalum is dispersive but 
does not show a distinct elastic-plastic structure. Although 
the final interface velocity is higher for the Lexan-backed 
impactor due to the incomplete unloading, the slope of the 
rarefaction waves are essentially the same. This result is 
significant since the attenuation results discussed below are 
for Lexan-backed impactors. 

It 

Lexan is a polycarbonate sheet made by General Electric Co. 
and has a density of 1.20 g/cc. An approximate equation of 
state determined for use in wave propagation calculations 
is c R a 31 u p + 15 u*. 
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Figure 16 Complete Wave Profiles, Tantalum 


If target thickness is large enough compared to .impactor thick¬ 
ness, the release wave will overtake the compressive wave and 
attenuate the peak stress, as shown in Figure 17. The initial 
stress was ^ 69 kbar and the impactors were backed with Lexan 
which gave only 'v* 80% unloading. The overtaking point can be 
calculated assuming ideal elastic-plastic compression and 
release waves, which gives: 




“"A 

C T””U 

_ h _ 



cj+u 

L p 


( 22 ) 


where is elastic release wave velocity 
(assumed to be 4.3 mm/us for 69 kbar). 


\ 
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< 

SECTION VI 

i 

I 

SPALL FRACTURE 

i 

i 

! , : 

Spall fracture by plate impact results from reflection of 

compressive waves from a relatively low impedance interface 

(normally a free surface) and subsequent wave interaction. 1 

Spall studies were carried out with both passive and active 

techniques. Passive methods involve the recovery and exam- 

' * 1 

ination of shock-loaded specimens. Metallographid examination 

* i 

establishes the type and degree of damage,'which can be corre¬ 
lated with test parameters such as velocity, impactor thickness 
and target thickness. Active methods utilize the laser velo- 

i , * 

city interferometer and provide time-resolved data;on the in¬ 
fluence of spall fractures on shock wave profiles, as measured i 
at the rear surface. . 


RECOVERY TESTS , ; 

i , 1 

Spall behavior of tantalum was studied by carrying out a series 

* ‘ ■ 

of impact and recovery tests, where the target was sectioned 
across a diameter, polished, etched, and examined optically a.t 
a magnification of 50X. The specimens were then graded or 
classified according to the degree of' fracture that was observed, 
which ranged from no visible damage to complete,material separa¬ 
tion. The onset of significant fracture is generally referred 

! 

* TJol 

Hanneman studied the thermal spall behavior of tantalum- 
using a capacitor discharge technique. He reported a "spall 
threshold temperature" of 1900 to -2140°C f6r tantalum wire. 

* " I 

) 

I 
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. to.as incipient spall and is a critical, point since it can be 
used to deduce dynamic fracture strength, and corresponds to 
generation of sufficient free surface area within the material 
ito reflept a portion of the interacting release waves as a 
.compiressive wave. The incipient spall threshold was defined 
as the impact velocity .(for a given set of impact parameters) 
corresponding to crafcking over at least 50% of the width of 
the section estimated to be under a condition of plane strain 
' during the time of loading* 

' • i i 

Incipient spall velocity results are summarized in Table IV 
'and the data are given in Figure 18. The symbols are defined 

i i ; 

ass i 

• , f • Complete Separation 

0 Above Incipient 

i . J ® Inqipient Spall '• 

9 iBelow Incipient 
0 No Visible Damage 



, Figure 18 Spall Data for Tantalum 
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As is typically the case with metals , the impact velocity for 

incipient spall increases for decreasing impactor thickness. 

This implies that peak stress (and therefore strain) required 

to create spall fractures increases with decreasing pulse 

width or time of loading. Optical photomicrographs of 

* 

recovered specimens are shown in Figures 19 and 20. The speci¬ 
mens were polished and etched to provide maximum contrast 
between cracks and sound material. The procedure followed 
was: 

1. Wet grind/ 320/400/600 grit silicon carbide. 

2. Rough polish, 6 pm diamond with 6 parts alcohol and 
1 part Gama-Medica Surgica Soap (Huntingdon Labs). 

3. Final polish/ 1 pm diamond with alcohol and soap. 

4. Swab etch (20°C), equal parts nitric acid, hydroflouric 
acid and water. 



0.185mm/ps 0.191mm/ps 0.213mm/ps 


Figure 19 Spall Fractures in Tantalum 
(1.5 mm 4.5 mm) 


All spall photomicrographs in this report are oriented such 
that initial shock wave propagation was from bottom to top. 
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TANTALUM SPALL RESULTS - 2'i°C 
0.41 mm-*-1.27 mm 1.50 mm-*-4.54 mm 





0. 213 mm / u s 

Figure 20 Spall Fractures in Tantalum 

As noted above, spall cannot be rigorously defined in terms 
of a unique impact velocity or stress, but requires considera¬ 
tion of the degree of fracture. Fracture in tantalum is of 
a ductile nature with the development of spherical voids. As 
velocity increases (see Figure 19), the voids begin to 
coalesce and a complete fracture surface is formed. Note that 
there is dispersion of the voids around a nominal "spall plane". 

This was also found for 6061-T6 aluminumand titanium, 

(52) 

i v* have well-defined elastic wave structures, while copper' 
has a poorly defined, low amplitude elastic wave and showed 
very little void dispersion around the spall plane. 

Scanning microfractographs of fracture surfaces for completely 
spalled specimens are given in Figures 21 and 22. Fracture 
was by normal rupture with evidence of extensive plastic flow. 

At -195°C, fracture is still ductile with no evidence of a 
ductile-brittle transition, but void size is much smaller. 
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Figure 21 Spall Fracture Surfaces in Tantalum 
(1.5 mm 4.5 mm, 0.32 mm/ys) 



lOOym 


Figrue 22 Spall Fracture Surfaces in Tantalum, 
Temperature Effects 
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SPALL PROFILES 


Complete wave profiles were obtained under spall-producing 
conditions using a tantalum-into-tantalum configurat-l'Mi, with 
a free rear surface on the impactor to give complete release 
and a free rear surface on the target for reflection of the 
compressive wave. The wave interactions resulting from such 
a test are discussed in Reference 53. Measured profiles for 
unattenuated pulses are given in Figure 23, with photomicro¬ 
graphs of the recovered and sectioned targets. Also shown for 
reference are the incipient and complete spall levels for this 
impactor/target combination as determined from recovery tests. 
Evidence of fracture is shown as reversal of the release wave. 

The amount of pullback (decrease in free surface velocity to 
the point of first reversal) may be related to spall strength 

of the material for a given impact geometry. An empirical 

(54) 

relation for this has been given by Taylor as: 

a s = pC(Au fs /2) (23) 

where o g is spall strength, p and C are local 
density and longitudinal wave speed, and Au^ 
is pullback. 

The data in Figure 23 give o g * 68 kbar. This calculated spall 
strength is not constant for a given material since it is 
proportional to Au^. This factor has been found to change 
with target temperature, pulse width and pulse shape for 
'1-T6 aluminum,and other metals such as tantalum may 
s..jw similar behavior. A spall profile for polycrystalline 
tantalum reported by Taylorshowed a pullback of 0.12 mm/ysec 
for a pulse width of ^ 1.5 ysec, which is less than the pull¬ 
back of 0.19 mm/ysec shown in Figure 23 for a 0.8 ysec pulse 
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width. However, for a given set of impact conditions 
(X Q , X, temperature), the pullback is relatively insensitive 
to maximum compressive stress as well as to degree of fracture. 
The spall profile from an attenuated pulse test, where the 
impact velocity was above that required for spall, is given 
in Figure 24. 


SPALL LEVELS 



Impact 


Complete 




0.8 1.2 



0.8 

TIME < sec 


Figure 23 


■ .'v, 

«••••'• ; ;• •* x —-> 

V-X 


! 
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Figure 24 


Attenuated Spall Profile i n Tantalum 
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SUMMARY 


The dynamic properties of tantalum were measured under uni¬ 
axial stress and uniaxial strain conditions. The material 
tested was commercially-pure bar stock with 45 to 500 urn grain 
size and 52 hardness. Material response to compressive uni¬ 
axial strain deformation approximated elastic-perfectly plastic 
behavior and exhibited strain-rate sensitivity with *v» 1.4 kbar 
yield at 0.001/sec and 'v* 5.4 kbar yield at 800/sec. 

Longitudinal and shear wave velocities were measured as func¬ 
tions of temperature and pressure, giving the following (in 
mm/visec): 

C r = 4.149 - 0.000146T 

J-» 

C s = 2.036 - 0.000173T 

(T in °C) 

C T = 4.146 + 0.00291P 

Jj 

C s = 2.032 + 0.00137P 

(P in kbar) 

Principal elastic constants w°re evaluated at T = 20°C and 
P = 0, including: 
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Bulk Modulus = K® = 1946 kbar 


! 

f 3K S \ 

[Wj 

S 

3.77, 

(*) 

T = 3.79 

Shear Modulus = 

G® = 688 

kbar 


1 

) 

ss 

s 

1.25, 

(#) 

T = 1.28 

Poisson's 

Ratio 

= 0.342 




Elastic wave velocity data were also used in determining the 
Gruneisen parameter y and the Debye temperature 0: 

y = 1.63 @ T = 20°C and P = 0 

= - 0.00043/°C 6 T = 20°C 

= - 0.00077/kbar § P = 0 
0 = 260°K 


The shock wave hugoniot was determined experimentally to be: 

a„ = 4.0 + 550 u +215 u* (a„ < 200 kbar) 

H p p H 

Assumption of an ideal elastic-plastic wave structure with a 
20 kbar hugoniot elastic limit and equilibrium initial and 
final states gave: 


49 








MANUFACTURING DEVELOPMENT 6 GENERAL MOTORS CORPORATION 


MSL-70-23, Vol. VI 

a H = 7.2 + 1820 m + 3120y 2 

P H = 1820 y + 3120 m 2 


U„ = 3.36 + 1.239 u 
S p 


The hydrostat as determined from elastic constants data was: 

(*)' 

(Keane equation of state) 

Measurements of compressive wave development showed the wave 
front to be characterized by a well-defined elastic portion 
and a transition to a spreading plastic wave. Elastic pre¬ 
cursor decay was 'u 50% to 19 kbar after 2 mm of travel. 

Spall fracture behavior was determined for elastic pulse widths 
of 0.20 and 0.72 usee. The incipient spall threshold expres¬ 
sed as impact velocity increased with decreasing pulse width. 
Tantalum showed a poorly-defined spall plane with appreciable 
dispersion of almost spherical voids. The fracture surface 
for completely spalled tantalum showed ductile failure with 
void growth and coalescence accompanied by a large amount of 
plastic flow. 


P„ = 6063 

n 


w 


1.10 


- 4723 In 


H 


+ 1261 
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APPENDIX A 

GRUNEISEN PARAMETER ESTIMATES 


1. The Gruneisen parameter is expressed thermodynamically as 



For an isotropic elastic solid, this becomes 



Substituting appropriate values for tantalum gives 


Y 0 = 1.63 

2. Slater developed a relation that was an extension of the 
Debye Theory for a three-dimensional continuum, and 
assumed Poisson's ratio was independent of volume and 
Y is the same for all vibrational frequencies: 



which simplifies at zero pressure to 
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where S is from U„ = C n + Su . 

on p 

Ultrasonic data gives S = 1.193 and y = 1.72. 
Shock wave data gives S = 1-239 and y = 1.81. 


3. Dugdale and MacDonald give an expression for plane wave 
velocity vs. volume in a three-dimensional continuum that 
leads to an expression for y differing from Slater's by 
a constant (1/3): 



which simplifies at **»ro pressure to 



From ultrasonic data, y = 1.39. 
From shock wave data, y = 1.48. 


4. 


Anderson and Dienes give expressions for y in terms of the 
longitudinal and shear modes of wave propagation: 
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Ultrasonic data gives y l = 1*69 and y s = 1.63. 

5. Schuele and Smith give expressions for the "low-tempera¬ 
ture" and "high-temperature" Gruneisen ratio, where low- 
temperature is the range in which it is assumed the solid 
behaves as a continuum and high-temperature is the range 
where all vibrational modes are excited and y is constant: 



Using ultrasonic data and y l and y s from above gives 
Y lt = 1.63 and y ht = 1.65. 


59 









MANUFACTURING DEVELOPMENT • GENERAL MOTORS CORPORATION 

»•* •« 

MSL-70-23, Vol. VI 

6. Perry calculated the Grueisen parameter from laser inter- 

£ 

ferometer data obtained for electron beam energy deposition, 
for polycrystalline tantalum, he obtained y = 0.95 assuming 
instantaneous deposition and y = 1.60 assuming 50 nsec deposi¬ 
tion time. 

(All ultrasonics data used in this Appendix are for zero pres¬ 
sure and 20°C). 


£ 

Perry, F. C., "Thermoelastic Response of Polycrystalline 
Metals to Relativistic Electron Beam Absorption", J. Appl. 
Phys. , Vol. 41, p. 5017-5022, 1970. 
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APPENDIX B 

DEVIATORIC STRESS ESTIMATES 


The shock wave equation of state for tantalum is given by: 

a = 4.0 + 550 u„ + 215 u 2 
H P P l 

a H = 7.2 + 1820 y + 3120 y 2 ! (B-l) 

P H = 1820 y + 3120 y 2 (Hydrostat) ' 

where y = (v q /v) - 1. 


The deviatoric stress (a D ) can be defined as the difference 
between the dynamic stress-compression path (a R - y) and the 
hydrostat, and is a measure of the material's ability to 
support shear. If the material strain hardens, the deviatoric 
stress will be a function of y. In Equation B-l, a D was 
assumed constant at 7.2 kbar which gives P„ = a„ - 7.2. 

n n 

The deviatoric stress can also be estimated from uniaxial 
stress data, shock wave data or ultrasonics data as follows: 



where Y is flow stress under uniaxial 
stress conditions. 
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A. The flow stress at low strain rate from the uniaxial 
stress tests (see Figure 3) is ^2 kbar, which gives 
a D = 1.3 kbar. 

B. For a hugoniot elastic limit (a g ) of 20 kbar, the rela¬ 
tion Y = a g (l-2v/l-v) gives Y = 9.6 kbar. From Equa¬ 
tion B-2, o D = 6.4 kbar. 

C. Since a can also be taken as the difference between 

cr„ and P„ at a compression corresponding to the hugoniot 
H xi 

elastic limit, one has 



For = 20 kbar: y e = 0.007 and P e = 13,5 kbar. This gives 
a D = 6.5 kbar. 

The estimates from B and C, above, agree reasonably well with 
the value of 7.2 kbar from Equation B-l, considering the assump¬ 
tions underlying each method and the experimental inaccuracies. 
The estimate from A, however, is significantly lower, which may 
be a consequence of the high strain rate sensitivity of tantalum. 
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Also, it should be noted that the stress-axis intercept is 
larger in the 0 H - y plane than the o H - Up plane. Assuming 
linear elastic and hydrostatic behavior up to the elastic 
limit, then: 



0 H - u plane - 

/ 2 2 
°D u “ rRT ^ P 0 C L ~ P 0 C B 

Using the approximations u g ^ u e /C L and P e - u e /(l+u ), 
and subtracting: 


4o D = °D 


“ C 0 C B U e (C L " C B ) 


For o e = 20 kbar, u e = 0.007 and A 0 Q =2.9 kbar, compared to 
Ao D =3.2 kbar from o H - u p and o H - y given in Equation B-l. 
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APPENDIX C 

EQUATIONS OF STATE, ULTRASONIC 


Numerous analytical and empirical relations have been developed 
for determining equations of state from ultrasonics data. Some 
of these are presenter here. 


1. The Birch equation of state for the isotherm can be 

* 

expressed as: 


3 T 

P zz ZL Tc 

T 2 0 




5/3 


where £ = 


12-3 K 


OT 





/ \ 7/3 / \ 5/3 


| 

[ / \ 2/3 


P T = 2891 

1 1+U ) - 1 1+P 1 


1-0.15 i 
j 

( l+y ) -l 



\ / ' ' 


L 1 

I v / 



Another parameter sometimes used with this equation of state 
is 



However, this term can usually be. neglected, contributing 
<0.4% at 100 kbar for tantalum. 
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2. The Murnaghan equation of state is derived by assuming 
the bulk modulus is linear with pressure which gives: 


K 


Isotherm - P = 


K 


0_ 
T' 
OT 


W 


K 


OT 


- 1 


P T = 507.1 


H 3 ' 


80 


- 1 


Adiabat - 


p 

S K S ' 

K os 


(■•■) 


K 


S' 

OS 


- 1 


P g = 516.2 


H - 


77 


- 1 


3. Instead of a linear bulk modulus-pressure relation, Keane 
assumed: 



S' S' 

where is the limiting value of K at very high pressure. 

This leads to: 
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(42i 

Anderson' ' suggested that this equation be applied by 

S 1 

empirically determining . This can be done by evaluat¬ 
ing the equation at a known high pressure shock wave point, 

S 1 

i.e., the equation is solved for K at a known P_, p 

00 O 

point. A high pressure data point was obtained by 
* 

Shipman at P fI = 1863 kbar and p = 0.497. This value 
of P H corresponds to P g = 1817 kbar, using the entropy 

correction term suggested by Duvall for the Murnaghan 

. S 1 

equation. Solving the above equation for K ro £t 

P g = 1817 kbar and p = 0.497 gives 



This givess 



Shipman, F. H., Unpublished Data, General Motors Corpora¬ 
tion, May, 1971. 

Huang, Y. K., "Analytical Approach to the Shock Compres¬ 
sibility of 18 Cubic-Lattice Metals", J. Chem. Phys., 

Vol. 53, p. 571-575, 1970. 
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APPENDIX *D 

WAVE PROFILES 


(Quartz Gage and Velocity Interferometer) 
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I 
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Tantalum, 5.0. 8mm Stock 
Velocity 3 0.261mm//tis 
Thickness= 6.11mm 
Tilt = 22ns 


0.6 

TIME (#isec) 


Tantalum, 76. 2mm Stock 
Velocity = 0. 263mm //js 
Thickness 3 6.30mm 
Tilt - 17.,s 


G.4 0.6 

TIME ifisez) 
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Tantalum, 50.8 mm Stock 
Velocity ■ 0. 362mm//is 
Thickness= 6. 10mm 
Tilt = 3Gns 


TIME (usee) 


Tantalum, 50. 8mm Siock* 
Velocity = 0. 258mm/jis 
Thickness = 4.58mm 
Tilt - 31ns 
Ihr. @ 1200°C 


TIME (usee) 


















INTERFACE VELC 





























INTERFACE VELOCITY (mm//isec) 
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X 

' "■ 

• • 

•' ■ 

‘ l 


1 

■ 1 
. „ . _i 

; : 



J 




TIME ( sec) 


0.6 

TIME (Msec) 


VI -154 

X -1.488 mm Ta 
o 

X - 4.552 mm Ta 
Vj - 0.270 mm(tsec 
FREE REAR SURFACE 


X Q -1.483 mm Ta 
X - 4.417 mm Ta 
Vj - 0.218 mmksec 
FREE REAR SURFACE 









FREE SURFACE VELOCITY (mm/psec) FREE SURFACE VELOCITY (mm//#sec) 
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